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Abstract 

We compare a momentum space implicit regularisation (IR) framework with other renor- 
malisation methods which may be applied to dimension specific theories, namely Differ- 
ential Renormalisation (DfR) and the BPHZ formalism. In particular, we define what 
is meant by minimal subtraction in IR in connection with DfR and dimensional renor- 
malisation (DR) . We illustrate with the calculation of the gluon self energy a procedure 
by which a constrained version of IR automatically ensures gauge invariance at one loop 
level and handles infrared divergences in a straightforward fashion. Moreover, using the 
tpf theory setting sun diagram as an example and comparing explicitly with the BPHZ 
framework, we show that IR directly displays the finite part of the amplitudes. We then 
construct a parametrization for the ambiguity in separating the infinite and finite parts 
whose parameter serves as renormalisation group scale for the Callan-Symanzik equation. 
Finally we argue that constrained IR, constrained DfR and dimensional reduction are 
equivalent within one loop order. 
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1 Introduction 



It is well known that the higher the symmetry degree of a quantum field theoretic model 
the more stringent are the constraints on a consistent regularisation scheme to handle the 
divergences which appear in diagrammatic expansions. For example, whereas a sharp cutoff 
may be successfully employed in most scalar field theories to reflect the correct physics in 
perturbation theory, it does not work so well already for abelian gauge field theories. For 
gauge field theories, DR is one of the most suitable frameworks because the amplitudes can 
be renormalised using for instance a minimal subtraction scheme (MS) and readily satisfy the 
Slavnov- Taylor identities. 

However some symmetries which are present in the integer dimension may not have a direct 
analogue in n dimensions. This is the case of supersymmetric (SUSY), chiral and topological 
field theories (the so called dimension specific theories). Some modifications of DR can be 
effected in order to mend certain shortcomings. For instance one may construct an extension 
of the algebra of the 7 s matrix to dimension n and control eventual spurious anomalies by 
imposing the Slavnov- Taylor identities as constraint equations. This is the usual procedure in 
the Electroweak Sector of the Standard Model. In Chern-Simons theories it may be necessary 
to employ an hybrid regularisation procedure by adding Higher Covariant Derivative terms 
in the Lagrangian which improves the ultraviolet behaviour. The remaining divergences are 
dealt with DR and adopting an extension of the Levi-Civitta tensor algebra to be compatible 
with analytical continuation on the space-time dimension. The main drawback in the example 
above is that the calculation may become extremely complicated beyond the one loop order. 
A variant of DR called Dimensional Reduction (DRed) was proposed by Siegel [|T| . The latter 
differs from DR in the sense that the continuation from 4 to n dimensions is made by compacti- 
fication. Thus whereas the momentum (or space-time) integrals are n- dimensional, the number 
of field components remains unchanged. Such procedure, however, may introduce ambiguities 
in the finite parts of the amplitudes as well as in the divergent parts in high order correc- 
tions. DRed has been largely employed especially in supersymmetric models as the invariance 
of the action with respect to SUSY transformations holds in general only for specific values 
of n. Unfortunately DRed appears to work well only at one loop level. In fact, DRed can 
be shown to be inconsistent in general with analytical continuation 0, || when 75 matrices 
and e Ml(U2 ... tensors are considered. In general a pragmatic attitude is adopted in handling the 
shortcomings brought by flawed regularisation frameworks especially when the model in con- 
sideration is known to be free of anomalies. In other words the task of treating the infinities 
in diagrammatic calculations, especially for theories which are sensitive to dimensional contin- 
uation, without introducing ambiguities steming from the regulator employed (that is to say, 
a regulator independent method) is still a subject of major interest. Ultimately it is desirable 
to construct a framework in which one has simultaneously : 1) no need to add structure to 
the Lagrangian and hence complicate the Feynman rules; 2) (nonabelian) gauge invariance is 
systematically guaranteed without having to be imposed as constraint equations order by or- 
der; 3) control upon infrared divergences without introducing additional machinery and 4) a 
method that is friendly from the calculational viewpoint. 

The task of treating the ultraviolet infinities in a regulator free fashion has been firstly con- 
ceived within the BPHZ formalism This framework relies ultimately in Weinberg's theorem 
which states that a Feynman graph converges if the degree of divergence of the graph and all its 
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subgraphs is negative. A systematic implementation of this idea is the Dyson's scheme which is 
based on the idea that differentiation with respect to the external momentum turns the graph 
less divergent. Hence in Dyson's method the divergent parts of a graph G are subtracted by 
applying Taylor operators t d ^ where d("y) is the degree of superficial divergence starting from 
the smallest subgraphs. When overlapping divergences occur care must be exercised in such 
subtraction procedure. The BPHZ framework is the generalisation of the Dyson procedure to 
include overlapping diagrams by means of a well prescribed formula called the forest formula. 
Although BPHZ is a very powerful framework which enables to construct proofs of renormaliz- 
ability to all orders, gauge invariance and hence the Slavnov Taylor identities should be imposed 
as constraint equations. The reason why gauge invariance is broken when the BPHZ method 
is applied to nonabelian gauge theories lies in the subtraction process which is based on ex- 
panding around an external momentum and thus modifying the structure of the corresponding 
amplitude. Some modifications in the BPHZ framework (Soft BPHZ Scheme) must be made to 
handle infrared divergencies because in the original formulation the subtraction is constructed 
at zero external momentum ||. 

DfR f7|-[]16|| and IR (please see []17 -p3~ for applications) seem to be very promising in this 
sense since they do not modify the space time dimension or introduce an explicit regulator at 
any step of the calculation. The former is position space method (contact with momentum 
space is made by means of Fourier transforms) whereas the latter is essentially constructed 
in momentum space. We shall discuss these methods in greater detail throughout this paper. 
We believe that the comparison which we shall outline here will show that IR is a promising 
candidate for handling divergences in field theoretical calculations (UV and infrared) in general 
in a symmetry preserving fashion yet being simple from the computational point of view. 

This paper is organized as follows: in section 1) we give a brief description of DfR and IR 
and compare with DR. We work out a few examples in 4 theory and QED where we discuss 
the role played by momentum routing invariance in connection with gauge invariance to effect 
a constrained version of IR. We also claim that to one loop order dimensional Reduction, DfR 
and IR are equivalent and define what is meant by minimal subtraction in IR. In section 2) we 
compute explicitly the setting sun diagram in both BPHZ, IR and compare with DfR. This is 
a nontrivial example because it possesses an interesting divergent structure from which we will 
clearly see the advantages of applying IR and DfR especially in obtaining the finite part. In 
section 3) we calculate the gluon self energy in QCD within IR to show that it can consistently 
handle the infrared divergences as well as readily display the finite part expressed by a class of 
well defined functions. We conclude by outlining a few applications in which IR could be useful 
and perhaps more advantageous. 



2 Relationship between Differential, Implicit and Di- 
mensional Renormalisation 

DfR was introduced by Freedman, Johnson and Latorre M as a method of regularization and 
renormalization in coordinate (Euclidian) space. The idea is that the product of propagators is 
not a distribution and so it has no Fourier transform. In DfR, renormalization is the procedure 
which extend products of distributions into distributions by substituting bad-behaved expres- 
sions by derivatives of well-behaved ones |J which are understood as distributions, that is to 



2 



say, the derivatives are meant to act on test functions. It automatically delivers finite Green's 
functions (which are identical to the bare ones for separate points but behave well enough at 
coincident points) without introducing an intermediate regulator or counterterm. For instance, 
suppose that we havejB^j/B^ an amplitude proportional to the product of massless propagators 

A (x) ' 



47T 2 X 2 



X 



Although ([!]) is a well defined distribution its square is not. According to DfR we search for G 
such that 

1 = UG{x 2 ) (2) 

which also guarantees manifest Euclidian invariance. In solving such differential equation we 
gain arbitrary scales among which M, which is introduced for dimensional reasons in 

^=-1^, (3) 

can be shown to play the role of scale variable in the (Callan-Symanzik) renormalisation group 
equation satisfied by the renormalised amplitude. The latter is constructed by substituting the 
l.h.s. of (OH with its r.h.s. where G given by tffi) , that is 



A* Or) 



Att 2 x 2 



r 2\2 



1 lnx 2 M 2 

4 x 2 



(4) 



R Htt 

Now G{x 2 ) does have a Fourier transform, namely ir 2 In (p 2 /M 2 ) /p 2 which enables us to write 
the Fourier transform of (^) as (in the Minkowski space) 

1 



k k 2 (k- P y 



R 



where 



(4tt) 2 ' 



(5) 



(6) 



J k = J G? 4 fc/(27r) 4 and M = 2M/^e, 1e being the Euler's constant. A comparison between DfR 
and DR's can be easily made. For the sake of clarity, we briefly outline it here for the case of 
massless theories following |J, [15 



Power law singularities of the type |x| n cannot have their degree of divergence decreased 
by using the identity 

r-|| r |-p+2 

\ x \-v = u l x l (7) 

1 1 (-p + 2)(n-p) v ' 

and setting p = n because of the pole l/{n — p). Alternatively we may try and regulate by 
dimensional continuation moving away from n by —re and thus using (R]) to get 
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\x 



n-2 



+ 0{e) . 



Thus in the dimensional approach the singularity x = is regulated by an infinite counterterm 
proportional to 8^ n '(x)/e. There is also a O(e ) piece in the term proportional to 8^ n '(x) (finite 
counterterm). If we subtract such counterterms we will be left with a term which is just 
the result obtained within DfR after identifying M with //. Alternatively we can use (|8[) to 
compute the regularised value of Aq(x). Given that the massless propagator in n dimensions 
is AqO) = -r(n/2 - l)|x| 2 - n /(4vr n / 2 ) and nA (x) = 5 (n \x) we have 



2e T\n/2-l) x4 _ 2n 

' A 9 r> 



4% n (4vr 
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+ 0(e). (9) 



Some comments are in order: The set of rules of DfR [ID| which fix local counterterms to 
establish equation (Q) is called constrained differential renormalization (CDfR). In particular, 
in CDfR one does not introduce arbitrary constants for singular behaviour worse than x~ 4 . 
CDfR can be shown to implement gauge invariance automatically at least to one loop order. 
From (P) it is clear that CDfR and DR's with a fixed initial condition given by (H) are equivalent 
under a MS scheme upon the identification M 2 = h 2 7t^eg 2 - As a matter of fact CDfR is identical 
to dimensional reduction to one loop order. In dimensional reduction the coefficients of the 
basic functions (finite, non-counterterm parts) are never projected into n dimensions because 
all the algebra is performed in physical dimension of the theory just as in CDfR. This is not 
the case of DR in which n can appear multiplying the basic functions which in turn produce 
different results from dimensional reduction. 

IR is a momentum space framework which somewhat resembles BPHZ in the sense that 
one algebraically manipulates the integrand of the amplitude in order to isolate the infinities. 
The idea is to isolate the divergences as basic divergent integrals(independent of the external 
momenta), e.g. , 

ho 9 (m 2 )= L_. (10) 

N 



9K } ~ } k {k 2 -m 2 ^ 

by using judiciously the identity 

1 ^ (-If {kf + 2k t ■ k) 



LV „ hi j ; (k 2 -m 2 ) j+1 

(k 2 - m?) N+1 \{k + hf - m 2 ]' 

where ki are the external momenta and iV is chosen so that the last term is finite under 
integration over k. Such basic divergent integrals which characterise the divergent structure 
of the underlying model need not be evaluated: they can be fully absorbed in the definition 
of the renormalization constants. We shall come back to this issue in connexion with what is 
meant by MS within IR and its relation to DR and DfR. An important ingredient of IR is that 
local arbitrary counterterms parametrized by (finite) differences of divergent integrals of the 
same superficial degree of divergence may be cast into a set of consistency relations |H|,p]||. 



They were shown to be connected to momentum routing invariance in the loop of a Feynman 
diagram. Should they vanish (as indeed they do in DR) then one would automatically have 
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momentum routing invariance and (abelian) gauge invariance. In other words by setting the 
consistency relations to zero (say, constrained IR (CIR)) one has the analogue to CDfR at one 
loop order. For n = 4 they read: 

T 2 = f — 2 [ kjxku 2^ 
^~] k k 2 -m 2 J k (k 2 -m 2 ) 2 ' 1 ' 



yO — I 9jiv a f k^k v 

v» - I. (k 2 -m 2 ) 2 J k (k 2 -m 2 ^ 3 ' ' 



np2 _ / 1 O / k^kykgk/3 

,va? ~ 9{^9 a p} I k2 _ m2 » J k (jfe2 _ m2)3 , 114 ) 



k 



-y-- ( ) / /* ^fl ^Ct fcf3 / \ 

^ = 9 ^ 9 <*ft I. (k 2 -m 2 ) 2 ~ 24 J k (k 2 -m 2 Y' (15) 



k 



etc., where g{^g af 3} stands for g^gap + g^gup+g^gua- Generically we may write = a^g^, 
etc. with oti arbitrary and finite. 

It is well known, however, that a shift in k is immaterial only if A s < 0, A s being the 
superficial degree of divergence, otherwise a surface term should be added. This is an indication 
that one should be careful in what concerns the momentum routing when divergences higher 
than logarithmic arise in Feynman diagram calculations. Perturbation theory makes a peculiar 
use of this feature for in some cases gauge invariance relies on adopting a special momentum 



routing ||24)| . A related issue is that whilst a shift in the integration variable is allowed within 
dimensional regularization, the algebraic properties of 75 clash with analytical continuation on 
the space-time dimension. In such cases, in IR we work with arbitrary values for the consistency 
relations until the end of the calculation so that physical conditions determine (or not) their 
value. For instance, a democratic display of the Adler-Bardeen-Bell-Jackiw triangle anomaly 



can only be achieved for arbitrary values of (|T3|) . 
2.1 Examples 

Here we illustrate the correspondence between the different regularisation frameworks in the 
context of a MS renormalisation scheme in massless 4 -theory and QED. In particular we study 
the Ward identity involving the QED vertex function whose finite parte is easily obtained within 
IR and analyse the role played by the consistency relations and arbitrary momentum routing. 
The 4-point function of the g/4!04-theory to one loop order T\(p)is proportional to (/x is an 
infrared cutoff) 



A ^ l h (k 2 - fj,) 2 [(k-p) 2 - fj?) (16) 



In IR we apply (jTT|) once to get 



p 2 z(z — 1) 



A = 3I w (/i 2 ) - b / dzlni K L + l), (17) 



V h 1 



with p 2 = s,t,u. In (17) we must separate the ultraviolet and infrared divergences (for the case 



[i — > 0) before proceeding to renormalisation. This can be easily accomplished by using the 
identity 

I log ^ 2 ) = I log (X 2 ) + bln^, (18) 
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which holds for arbitrary A. This enables us to write 



rU?) = f(wA 2 ) + Hn(^)- 



1 

b j dzhi 
'q 



p 2 z(z — l) 

-V 2 



+ 1 



K _ 2 + Q> (^±i)), 



(19) 



with a s 
to yield 



v/v7 



■s + 1. Hence we define the MS within the IR method by subtracting Iiog{^ 2 
b g 2 f ^ f\ 2 \ , /a B + 1 

~~ 2 



1 R 



MS 



(P) 



In 



— a. In 



- 1 

M 2 



(20) 



2 \ v/i^/ \a, 

which is just the result obtained in DfR JTJJ with A 2 = A 2 e 2 = M 2 . Notice that in the limit 
where \i — > 0, /i cancels out in the equation above, as it should . That is because the infrared 

divergent piece of the logarithm b In Q?) cancels out with another piece coming from the finite 

part of the amplitude, this enables us to write 



A 



R 

MS 



-A 2 



(21) 



Therefore we have in the MS scheme as defined above for IR the same prescription as defined 
by equation (|5|) in DfR. Moreover one can verify that A plays the role of renormalization group 
scale in the Callan-Symanzik equation. This is expected since it parametrises the arbitrariness 
in separating the divergent from the finite part. In the massless limit the Callan-Symanzik 
equation for the 4-point function reads 



x 9 a 9 , 



o 



(22) 



from which we compute the standard value f3 = 3g 2 /(167r 2 ). 

We also expect IR to be identical to dimensional reduction (as DfR is) to one loop level 
for the Lorentz algebra which determines the coefficients of the finite parts in IR is effected 
in the integer dimension, say n = 4. In order to illustrate this point, as well as to pinpoint 
the role played by the consistency relations in IR in connection with momentum routing and 
gauge invariance let us study the QED Ward identity involving the vertex function in IR. The 
electron self-energy in the Feynman gauge is written as (e 2 = 1) 



1 



7m" 



k Jk+ Jkx — m (k + k 2 ) 2 - Ai 2 
| = (A-2m)/ + 7 ^ 



(23) 
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where /i 7 is an infrared regulator, ki,k 2 are arbitrary momenta running in the loop such that 
k\ — k 2 = p, p being the external momentum, which we shall parametrise by a by setting 
ki = (1 + a)p and k 2 = ocp), and 



7, I* 



hk>* 



Now within the spirit of IR we use ( |i~H ) to write 

f k\ + Ik ■ k 2 + m 2 - /I 2 

I = Ii g{m ) - 



k (k 2 - m 2 ) 2 [(k + k 2 ) 2 - fi: 
kf + 2k- ki 



k (k 2 + m 2 ) [(k + k 2 ) 2 - fi 2 ] [(k + h) 2 - m 2 } 

2 

7' 



I log {m 2 ) - b Z (jl, m\p 2 - m 2 ) , (24) 



in which Z is part of a class of functions which characterise Feynman diagram calculations to 
one loop order, 



Zk{^l, A 2 ,,]? 2 ; A 2 

;(l-z) + 

-A 



'^' l ./" 1 -" + 'l^-- A - (25) 



We can similarly calculate J M to get 



^ = + k 2 fI log {m 2 ) + (h + k 2 ) Aip" 

+ bk%Z {fi 2 m 2 ,p 2 -,m 2 )-bp»Z 1 { f i 2 m 2 ,p 2 -,m 2 ). (26) 



where Ai is defined from (|T3|) as 



T°uu = Ai57«/ 



and is in principle undetermined. The limit /i 7 — > is well defined for the Z k functions and 
Zi(0, m 2 ,p 2 ; m 2 ) = 1/2Z (0, m 2 ,p 2 ; m 2 ). This enables us to write 

S(p) = -(i>-4m)(/ w (m 2 ) -6Z (0,m 2 ,p 2 ;m 2 )) 

+ (2a + l)Ax £). (27) 

In order to establish the value of Ai it is natural to check whether the Ward identity which 
relates E to the vertex function can place any constraint on Ai. Consider the QED vertex 
function with incoming momenta p and q and outgoing momentum p + q Q: 

- ' A " M -b^^jrj^^j^, (28) 



1 Because the QED vertex function is superficially logarithmically divergent, it must be momentum routing 
independent. 
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Within the framework of IR we can write, after some tedious yet straightforward algebra, 

-iA"(p, q) = 7 ,T^ + ^I log (m 2 ) + 2bf{Z - Z ) + 
+ b£°° Um(p + g) M - 2mV - 2 ^ A + 

+ 26( J 6+^) 7 ^e 01 + 2M7 M (i'+^^ - 

- 8mb(jf £ 01 + g" £ 10 ) - 46(p" ^ 02 + <f M™) + 

+ Kp^ + ^K 11 , (29) 



where 



Z = Z (m 2 , m 2 , (p - q) 2 ; m 2 ) , 



Z = Z(fi 2 m 2 ,p,q) 



7' 

l rl-z /n(„ „ „, „ ,,2 ^ 1 2^ 



dz / c?y In 



r n = r n (^,m 2 ,p, q ) = 

1 /-1-2 y n z m 

dz / ' 



-m 2 



o ""Jo Q(P,Q,y,z,fi 2 ,m 2 ) ' 

<3(p, y, /4> ™ 2 ) = p 2 y( 1 -y) + i 2 ^ 1 - *) - 

-2p-qyz + (z + y)(fi 2 -m 2 )-fi 2 . (30) 

Now we write Tq U = X^g^ in (p9j), A2 being an arbitrary parameter, and redefine 

E(/) = S( i 6) + (2a + l)A 1 
A"(p,g) = A"(p,g) + 7*A a (31) 

which with the help of the relations displayed in the appendix, enable us to verify promptly 
that 

{p-q)^(p,q)=Z(ft)-Z{jq). (32) 
Hence the Ward identity is fulfilled if 

A 2 = (2a + l)Ai. 

The natural choice is to set Ai = A2 = 0, which automatically implements both gauge and 
momentum routing invariance (CIR). Notice that by setting Ai = in fl2"T|) leads to the same 
result as obtained in CDfR and dimensional reduction Jl6| (which is however different from DR 
within the same subtraction scheme) as we too have worked in four dimensions. 

This illustrates the equivalence between CIR, CDfR and dimensional reduction to one loop 
level. For instance, the superfield calculation of the one loop correction to the vector propagator 
is gauge invariant in dimensional reduction scheme only if |2(J 
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We can easily evaluate the integral above within IR to show that it reduces to 

= P"K» ( 34 ) 

showing that CIR (T° = 0) may be safe framework to handle the problem. In principle IR 
is generalisable to higher loop calculations avoiding breakdown of symmetries such as gauge 
invariance and supersymmetry [ |25|| . 

3 IR, BPHZ and DfR: A Two Loop Example 

In order to illustrate the correspondence between the BPHZ formalism and IR, we shall compute 
the 2-loop correction to the 2-point function in <fi\ theory in both methods and compare with 
DfR. The amplitude is depicted in figure (|l|) (setting-sun diagram). The computation of the 
finite part of this diagram is notoriously difficult in DR, for instance. However for both IR 
and DfR it can be readily displayed. The BPHZ scheme relies on the forest formula to 



k i + p 




-*1 - *2 



Figure 1: Sunset Diagram 



perform the subtraction of the divergences from an amplitude Let Iq be the integrand of 
such amplitude associated with a graph G. Then the subtracted integrand is given by 

Rg = J2U -< Ml G (35) 

where is the set of all the forests U of G, including the empty set Q, d{^) is the super- 
ficial degree of divergence of the subgraph 7 and t d ^ is the Taylor operator which cor- 
responds to an expansion around to order ^(7) in the external momentum to the sub- 
graph. For the sunset diagram the subgraphs are shown in figure Therefore we can 
write 4> = {0, G, 71, 72, 73, G71, G72, Gj 3 }. Thus 

R G = (1 - - t° y2 - t° 3 - t G + 

+_ 4< + tlA, + ^A^g (36) 

2 A forest U of a diagrama G is the set of all subgraphs of G, including G itself, which are neither overlapping 
nor trivial. 
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Tj ( k i = ) 



Figure 2: Subgraphs 71 , 72 and 73 



which are ordered so that if 71 C 72 then t 71 lies on the right of t 12 . The amplitude for the 
sunset diagram is superficially quadratically divergent. It reads 

JOO 

I™ = A(ki + p)A(k 2 )A(k l + k 2 ) with 
A W S (37) 



Therefore we get, using (|36|). 

R G = Ig - tll% - A{k 2 )\l - t 2 G ) A(* a + p) . (38) 

If one uses a regulator, it can be shown that the term A{k 2 ) 2 {l — t G )A{k\ +p) above actually 
vanishes. This comes from the fact that 



F(p)= f(l-t 2 G )A(k+p) 
Jk 



is independent of p and therefore F(p) = F(0) = 0. Thus according to the BPHZ forest formula, 
the finite amplitude is obtained through the subtractions: 

R ( ; = I I%- [ t%ir 



G— G / L G 2 G > 

fci,fc 2 Jk\,k.2 J ki,k 2 



I t 2 G Ig = I U l(fl 2 )-p 2 hl2^ 2 ) + 
J k\,k2 

+4 / (p • A; 1 ) 2 A(A; 1 ) 3 A(A; 2 )A(A; 1 + fc 2 ) , (39) 

J kl,k2 

where I mnp (/i 2 ) stands for 

(40) 



k lM (*? - fi 2 ) m (kj - ii^dh + k 2 y - ■ 
10 



The integrals I mnp are convergent only if m + p > 2, n + p > 2, m + n > 2, m + n + p > 4 are 
satisfied []. In order to display the finite part of the setting sun amplitude there is still some 
work to do in (|39|). In particular one would have to adopt a regulator to proceed in this task, 
the result being regulator independent as guaranteed by the BPHZ scheme. That is to say, the 
BPHZ scheme guarantees to us is which subtractions one ought to perform in order to render 
the amplitude finite and that the result is regularization independent. 

Now let us evaluate fl37|) in the light of IR. Because d{G) = 2 all we need is to apply (|TT| ) 
up to N = 2 in A(ki +p) so to display the infinities in terms of basic divergent integrals which 
are cleared out of external momenta to get 



A{h+p) = A(k 1 )-(p 2 + 2k 1 -p)A 2 (k 1 ) 
+ (p 2 + 2p-k 1 ) 2 A 3 (k 1 ) 
- {p 2 + 2p-k l fA\k l )A{k l +p) 



(41) 



which enables us to write 



hll(^)~p 2 hl2(^) 

+ 4 / (p • fc 1 ) 2 A(A; 1 ) 3 A(A; 2 ) A{h + k 2 ) + p 4 I 3U 
! (p 2 + 2p-k 1 ) 3 A(k 1 ) 3 A(k 2 )x 

J k\ ,k2 

x A(/c 1 + A; 2 )A(A; 1 +p). 



(42) 



Notice that the first three terms in the rhs of the equation above are just the terms which were 
prescribed to be subtracted using the BPHZ scheme; the fourth term and fifth term (let us call 
them F4 and F§ for definiteness ) are clearly divergent but their difference is finite. To see that 
let us isolate the divergence in both terms as a function of one loop momentum variable only 



using (|lTl) . Thus 



whereas 



F, = p 4 I loa ^ 2 ) [ A(ki) 3 - p 4 x 

Jk! 

x I (k\ + 2ki ■ k 2 )A(ki) 3 A(k 2 ) 2 A(k 1 + k 2 ) , 

Jki,k2 

F 5 = Ii og (v 2 ) [ (p 2 + 2p-k 1 ) 3 A(k 1 ) 3 A(k 1 +p) 

Jk! 

~ j (p 2 + 2ki ■ p) 3 (k 2 + 2ki ■ k^Aih) 3 x 

J ki,k 2 

xAihfAih + k^Aih+p) . 



(43) 



(44) 



3 It will also be useful to know that I!^L = f, , m — 2 \l?n \ ^^li^A'^l ^2 — jtz converges if 2m + 2p > 

mnp Jki,k 2 (k( — u 2 ) m (fc| — u z )™ ((fcx +k 2 ) z — ^ ) p ° 1 

4 + 2i + k,2n + 2p>4: + 2j + k 7 ni + n+ p>i + j + k + 4: and m + n>i+j + k + 2 are satisfied. 
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By using Feynman parameters to evaluate the first terms in F4 and F5 one can show that they 
cancel out. Therefore we can write the finite part T = F 4 — F 5 of the setting-sun amplitude as 



T = 7T / (V + 2A;i • P) 3 (*i + 2A;i • *&) x 
6 J kite V 

xA(A; 1 ) 3 A(A; 2 ) 2 A(A; 1 + k 2 )A(h + p) - 
-p 4 (A; 2 + 2An ■ fc 2 ) A(A; 1 ) 3 A(A; 2 ) 2 A(A; 1 + k 2 )) . 



(45) 

In order to make contact with other regularisation/renormalisation frameworks let us take a 
closer look at the divergent structure of T^ 2 \p), namely 

^(/ m (/x 2 )-p 2 / 112 (/i 2 ) + 
+ 4 / (p • k 1 ) 2 A(k 1 ) 3 A(k 2 )A(k 1 + k 2 )) . 

Jkite 

It is easy to show that the third term above may be written as 

P 2 (h^ifJ 2 ) + fJ 2 ImifJ 2 )) and that 
hia = -—IiogifJ 2 ) + N ,N finite 



(46) 



2^ 



N = - j (k\ + 2fci • k 2 ) 

J kite 

xA(A; 1 ) 3 A(A; 2 ) 2 A(A; 1 + A; 2 ) . 



(47) 



Also let us split the logarithmic divergence using Collecting all the results together we 

have 4 



<2)R 



b 2 P \ /A 2 
In — 



(p) = £ 
' 6 V 2 

(P 4 + + 2fci • A; 2 )A(A; 1 ) 3 A(A ;2 ) 2 x 



kite 



xA(h + k 2 )+ j {p 2 + 2/n • p) 3 (A; 2 + 2A* • fc 2 ) x 

J kite 

xA{k 1 fA{k 2 fA{k 1 + k 2 )A{k 1 +p)\ , 



where 



The equation ([ 
in DfR, 



rg?(p)-rg> fl (p) = ^(W)-^WA 2 )). 



(48) 



(49) 



obtained within IR is the momentum space analogue of the result presented 



111 



n? 



{x,M) 



96(4vr 



2^3 



:n-9// 2 )(n-// 2 ) 



J In CDfR the dimensionful constants are taken from the logarithmic divergent pieces only (lz 
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M 



+27r 2 ln^-(n-/i 2 )(5( 4 )( 



x 



(50) 



Here too, as it was claimed in Ol for DfR, one has been able to display the finite part of the 
setting sun diagram in a closed and compact form in an easy fashion with the advantage of 
working directly in the momentum space. Notice that the scale A in (48) is just the analogue of 
M in the equation above and plays the role of scale in the Callan-Symanzik equation satisfied 
by rjjj (p). For instance, for fi — > 0, T^Qo 2 ) is well defined and obbeys 



d 



d 







from which we may calculate the lowest order value of 7^, namely 



7<t> 



12 (167T 



2\2 ' 



(51) 



(52) 



In fact it can be shown that 7^ to lowest order is simply the coefficient of the logarithmic 
divergence lnA 2 //z 2 p7| . By using a general parametrisation for Ii og (X 2 ) p0| , viz. f] 



Ilog{>? 



bin 



+ ?7 



(53) 



(A is an UV cutoff and 77 an arbitrary constant) in ([4~9|) we see that the coefficient of the 
logarithmic divergence evaluates to 7^ given in fl5"2"|). Alternatively one can apply DR to evaluate 
I log (X 2 ) which gives &r(e) = 6(l/e - 7 £ + 0(e)). 

We can also study the setting sun diagram with arbitrary routing . Let us split the external 
momentum p between the upper and lower lines in figure (|J) so that the amplitude reads 

,2 



6 



A(fci + ap)A(k l + k 2 )A(k 2 + (3p) 



ki,k2 



9 -7T \ Ig{<*,0) 
Jkito 



(54) 



with a + (3 — 1. Using the forest formula one can show within the BPHZ framework that 



R, 



G 



I£(a,P) - I m (//) + (a + /3) 2 / U2 (/i 2 ) 



-4(a + /3) 2 I (p ■ fc^A^Afo) Afo + k 2 



(55) 



To work out ([54]) within IR, it suffices to expand A(ki + ap) and A{k 2 + f3p) using (jTTf) up to 
= 2 to get 

rg } (p,a,/3) = J m (/i 2 ) - (a + /3)V/m(/i 2 ) 



+4(a + /5) 2 

'fcl,/C2 

+ finite (a, (3 = 1 



(p-fci) 2 A(fci) 3 A(fc 2 )A(fci + fc 2 



(56) 



3 Such parametrisation is constructed based on dIi og (X 2 )/d\ 2 = —b/A 2 . 
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which turns out to be identical to (J£2D, as it should, because the finite part above can be 
shown to be independent of a. No consistency relations have appeared in this example as the 
momentum routing independence in the setting sun diagram is trivial. In [21] we calculate the 
j3 function of <fi\ theory to two loop order within IR. 



4 Gluon Self-Energy of QCD 

In both abelian and nonabelian gauge field theory the gauge boson self energy is bound by gauge 
invariance to have the structure i(p 2 g^ — P^Pu)^(p 2 )- For QCD, the cancellations that lead to 
this structure are more complex than for QED and relies on a gauge invariant regularisation 
framework to handle both the ultraviolet and infrared divergences. It is well known that 
adding a small mass for the gluon breaks gauge invariance although it may be safely done for 
the photon. 

It will be interesting to see how gauge invariance is implemented in IR for the gluon self 
energy in connection with the consistency conditions ( |T2"D -fli"5]). Let p be the external momen- 
tum. The diagrams that contribute to the gluon self energy to one loop order are well known: 
1) the gluon tadpole nj(l), 2)the gluon loop 11^(2), 3) the ghost loop 14^(3) and 4) the 
fermion loop 11^,(4). Hence the gluon self energy is given by 



K = n£(i) + n£(2) + n£(3) + n£(4) (57) 



The Feynman rules in momentum space for QCD can be found in any textbook p7j- ft is 
straightforward to show that 



n£(l) = -g 2 C 2 (G)5 ab 3 [ ■ 

Jk 



k 2 — /j 2 

= -3g 2 gflu C 2 (G)5 ab I quad (fi 2 ) (58) 
The gluon loop is equally simple to be displayed within IR. It reads 

K P) - \ S/^^ v^Tk^? • (59) 



where 



Using that 
we may write 



N n» = [ g w(p- kf + g P a (2k + pY + g^(-k-2p) p ] 
x [Sp(k - p) a + 9pa (-2k - pf + %(k + 2p) p ] 
= 2p il p v - + Puk^) - IQk^ky 

- g^[{p-k) 2 + {k + 2p)% (60) 
(p - k) 2 + (k + 2p) 2 = (k + p) 2 + k 2 + Ap 2 , (61) 



nj(2) = - l -g 2 C 2 {G)5 ab [{2p,p v -Ap 2 g, v )J{p 2 ^ 2 ) 

- 10(^J>V) + J^(pV))], (62) 
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whereas the ghost loop is given by 



nj(3) 



2 rdac rcbd 



k 2 — /i 2 [{k + p) 2 — jl 2 



-g z r c f 



X 



x(p + k)^ v = -g 2 S ab C 2 (GKp u J^p 2 ,fi 2 ) 



in which we have used the notation 



J, 



+ b 



flis 



6 L3 



\-- 2 {p 2 -» 2 )z 

6 p z 

+ ^(-p 2 + 4m 2 )Z 



J, 



-2p a e^ + -p,bZ 



with the Z n functions defined as in 



and 



Z Q = Z (fi 2 ,fi 2 ,p 2 ;fi 2 ^ 



Q (o) 



0(2) 



e (o) 



k^kjy 



(k 2 - /i 2 ) 3 



kuk v 



(k 2 -n 2 ) 2 
kjj,k v k a kp 

{k 2 



2U 



Hence we can write ([U|) as T° = g^hogil^ 



46^, etc. Then it follows that 



(63) 

(64) 
(65) 

(66) 

(67) 

(68) 
(69) 

(70) 



(71) 



n£(i)+n£(2) + nj 



X 



g 2 C 2 (G)S ab (p% v - Pf , Pv ) 
- &| + | (/ /o9 (/i 2 ) - bZ 



(72) 



The fermion loop contribution to the gluon self energy is identical to the vacuum polarization 
tensor of QED except for the colour and number of fermions (Nf) factors. It has been computed 
within IR elsewhere [|2U|| . Here we only write the result in which we have already subtracted 
Iio g (^ 2 ) (that is to say we have employed the minimal subtraction in IR) and set the consistent 
relations to zero (CIR): 



UfU) = —a 2 ^- 

^ ' 144vr 2y 2 



$ ab {PnPv - p 2 g^u) x 



(l2Z (m 2 f ,m 2 f ,p 2 ;X 2 )+4 



(73) 
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Now the limit where /i — > can be safely taken because using (|TBD and that 



one can show that 



and 



Z (M 2 ,/iV;/i 2 ) = ln^ + Z (/i 2 ,// 2 ,p 2 ; A 2 ) , (74) 



-&Z (// 2 ,/iV;A 2 ), (75) 



2 



Z (0,0y ; A 2 ) = ln^-2. (76) 
Let us also take the limit of massless fermions. Thus we have, 



n£(4) = j^- 2 g 2 ^ ab (P»P» - pV)[i21u(a 2 / p 2 ) - 4] . (77) 



Bringing all the results together enables us to write the complete gluon self energy to one loop 
order, after minimally subtracting in the sense of IR the divergence expressed by Ii og (\ 2 ) and 
setting Y^, to zero as 

<^ = v w - P V) >< 

((15C 2 (i2) - 6N f )\n(\ 2 /p 2 ) - 2C 2 {R) + 2^) (78) 
which is just the result obtained in DfR |12| after identifying A with the differential renormal- 



ization arbitrary scale M . Notice that both the massive and massless cases can be straightfor- 
wardly handled in IR. 

5 Conclusions 

In this paper we compared three frameworks namely differential (DfR), implicit (IR) and BPHZ 
regularisation/renormalisation which being strictly defined in the physical dimension of the 
underlying theory may overcome the problems which arise when applying dimensional regulari- 
sation and variants to dimension specific theories such as supersymmetric, topological or chiral 
quantum field theories. 

The purpose was to motivate the answer to a few questions related to the consistency and 
applicability of IR in handling infinities in Feynman diagram calculations, namely: 1) study 
how infrared divergences are treated within this scheme; 2) understand how (nonabelian) gauge 
invariance can be automatically implemented within a constrained version of IR; 3) define what 
a minimal subtraction renormalisation is in IR in analogy with dimensional and differential 
renormalization; 4) argue on the equivalence between IR, DfR and DRed to one loop order 
and 4) motivate IR as a practical and consistent tool for loop calculations in dimension specific 
theories. 

Since in implicit regularization the divergences are displayed in the form of basic divergent 
integrals it is natural to ask what is meant by minimal subtraction in such scheme. We have 
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shown that the logarithmic divergences expressed by Ii og (n 2 ) can be split as in equation 
to give rise to am arbitrary scale which plays the role of renormalisation scale in the Callan- 
Symanzik equation. By subtracting Ii og (X 2 ) when a logarithmic divergence occurs we have a 
finite part which is identical to the result in differential renormalization (with A playing the role 
of the arbitrary scale M in DfR) and dimensional regularisation (except for a local counterterm) . 
We define such renormalisation scheme in IR as minimal. In constrained the arbitrary scale is 
also taken from the logarithmic divergences only ]15[ (functions with singular behaviour worse 
than logarithmic (a; -4 ) are reduced to derivatives of less singular functions without introducing 
extra dimensionful constants. This is the main ingredient that fixes the renormalization scheme 
in (constrained) DfR and automatically preserves gauge invariance at least to one loop order. 
We showed in the calculation of the gluon selfenergy that a constrained version of IR preserves 



gauge invariance just as it does for abelian gauge theories |19| . Constraining IR amounts 



to set some well defined finite differences between divergent integrals of the same degree of 



divergence to zero [19|], @- Such differences are called consistency relations and were shown 
to be connected to momentum routing invariance in the loops. 

In order to illustrate the relationship between the BPHZ and the IR schemes we have 
computed the setting sun diagram. The BPHZ scheme is a very powerful tool for all order 
proofs, for instance, and it delivers unambiguously the terms which ought to be subtracted in 
order to render the amplitude finite by means of the forest formula. It is a subtraction method 
which is regularisation independent in what concerns the finite part. In order to proceed the 
calculation so to extract the finite part, one has to apply ultimately a regularization scheme. 
However symmetries may be broken in the course of such operations such as gauge invariance. 
As we have seen within IR certain surface terms are important in order to preserve gauge 
invariance. Moreover, an expansion around zero in the external momentum potentially breaks 
the gauge invariant structure of the underlying amplitude. By controlling surface terms and 
using an identity at the level of the integrand flTTp we have verified that IR has control upon 
gauge invariance at least to one loop order. In other words in IR, the finite part is delivered 
automatically and no damage is made to the integrand whereas arbitrary local terms are duly 
parametrised in IR by the consistency relations. A proof of renormalisability to all orders in 
an alternative fashion to BPHZ has been constructed for tp\ theory within IR. 

Although we have restricted ourselves to the gluon selfenergy, surely one should verify if the 
other Slavnov- Taylor identites are satisfied as well, by calculating the vertex functions of QCD 

El- 

Finally we conclude that both DfR and IR are potentially good frameworks to apply in 
dimension specific problems in order o avoid ambiguities and spurious anomalies. IR is par- 
ticularly friendly from the calculational viewpoint with the advantage of working directly in 
the momentum space. We believe that computations beyond one loop order in Chern-Simons 
matter theories as well as in supersymmetric models may profit from IR since such method does 
not modify the underlying theory and operates in the physical dimension of the the theory. 



Appendix 

The following relations between the functions Zk and £ mn can be easily checked and simplify 
the explicit verification of the Ward identity involving the electron self energy and the vertex 
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function. 



? 2 6i -P-Qdo = ^{Zo(q 2 ;m 2 ) - Z (p.q;m 2 ) + p 2 £oo} 



(79) 



2 lf -Z (p + g) 2 ;m 2 ) Z (p 2 ;m 2 ) 
= 2 g + ^ 



+ 9^ 



10 



9 6)2 - = 

p 2 6o - p-g^n = 



+ m 2 £, 



oo 



+ m 2 £ ( 



oo 



+ ±Z (q 2 ;m 2 )+p% i y 



(80) 

(81) 
(82) 



(83) 



2 C 00\ 



Z(/i 2 , m 2 ,p, q) = l/2Z„((p - g) 2 , m 2 ) - (1/2 + ^ 
+ l/2(g 2 + /i 2 - m 2 )e° + l/2(p 2 + /i 2 - m 2 )£ 01 . (84) 

where £ mn = £ mn (m 2 ,m 2 ,p,q) and we have abbreviated 

Z k (m 2 ,m 2 ,p 2 ;m 2 ) = Z k (p 2 ;m 2 ). (85) 
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